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Highly-damped quasi-normal frequencies for piecewise Eckart potentials 

Jozef SkakalcE and Matt VisseiQ 
School of Mathematics, Statistics, and Operations Research; Victoria University of Wellington, Wellington, New Zealand. 

(Dated: 4 May 2010; 15 June 2010; IMfeX-ed September 2, 2010) 

Highly-damped quasi- normal frequencies are very often of the form uj n = (offset) + in (gap). 
We investigate the genericity of this phenomenon by considering a model potential that is piece- 
wise Eckart (piecewise Poschl-Teller), and developing an analytic "quantization condition" for the 
highly-damped quasi-normal frequencies. We find that this cj n = (offset) + in (gap) behaviour is 
generic but not universal, with the controlling feature being whether or not the ratio of the rates 
of exponential falloff in the two asymptotic directions is a rational number. These observations are 
of direct relevance to any physical situation where highly-damped quasi-normal modes (damped 
modes) are important — in particular (but not limited to) to black hole physics, both theoretical 
and observational. 
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In many diverse branches of physics one is interested 
in studying potentials that enjoy suitable falloff con- 
ditions at spatial infinity, and it is generally observed 
that such potentials lead to quasi-normal modes (QNMs, 
damped modes) with associated quasi-normal frequen- 
cies (QNFs). To help gain semi-analytic understanding 
of this phenomenon we investigate the QNFs of a piece- 
wise Eckart (Poschl-Teller) potential [iHj]. We are par- 
ticularly interested in understanding the 



tj n = (offset) + in (gap) 



(1) 



behaviour that has been encountered in very many differ- 
ent analyses, often in the context of black hole physics, 
but by no means limited to black hole physics [5r|38(. 
The specific model we are interested in is 



t/)"(x) + V(x) 4>{x) = 0, 



with a piecewise Eckart potential [3£ 



V(x) 



V - 
Vb- 



sech (x/b-) 
sech (x/b+) 



for x < 0; 
for x > 0. 



(2) 



(3) 



The standard case is Vq- = Vb+ = Vq, with fo_ = b + = b, 
so that V(x) = Vq sech (x/b). A related model where 
Vo_ = Vb+ = Vo but with b + ^ &_ has been explored 
by Suneeta [40|, but our current model is more general, 
and we will take the analysis much further. We start 
by imposing quasi-normal boundary conditions (outgoing 



radiation boundary conditions) [39, 41 



ip+(x — > +00) 



yb_(x — > —00) — > e 



+ IUJX 



(4) 



On each half line (x < 0, and x > 0) the exact wavefunc- 
tion (see especially page 405 of |4lJ) is: 



n (•'■) =' = ''"-'-' -iFi ( ^ + a±,--a ± ,l + ib±uj,z 



(5) 
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where 



\ - V b 2 for V b 2 < 1/4; 
V b 2 -\ for V a b 2 >l/A; 



(6) 



and z = 1/(1 + e ±2x / b± ). The key step in matching 
these two exact wavefunctions at x = is to calculate 
the logarithmic derivative. Using the Leibnitz rule and 
the chain rule one evaluates i/j' ± (0)/ijj±(0) as: 



=Rw=F 



1 din {2-F1 (i + a±, \ — a±, 1 + ib±uj,z)} 



2b4 



Az 



2=1/2. 

(7) 



Invoking the differential identity 
d{ 2 -Fi {a,b,c,z)} 



dz 
c-1 

z 



[ 2 Fi (a, b, c - 1, z) - 2 Fi (a, b, c, z)} , (8) 



=pza; 



1F1 (| + a±,i -a±,ib ± L 



we see 

? — 7 = ~\~ LUJ " / 1 1 \ 

V±(0) 2F1 (5 +a±, %-a±,l + ib±u,z) 

Now using Bailey's theorem 



F\ a, 1 — a, c, 



1 



r(§)iW) 



c+a \rf c-a+l \ ' 



we have the exact result 



V4(Q) 



T 



2 r( 



ct±-\-iusb± 



r(s±a)r( 

I)r( 



2=1/2. 

(9) 



(10) 



3\TV -a± +i<^h± , 3\ 
2 ^4^ 



5± r( "±+^± + l)r( -"±+^ b ± + I)' 



(11) 



To obtain a more tractable result it is extremely useful 
to use the reflection formula 

7T 



r(z)r(i-z) = - 



sin(7rz) ' 



(12) 



derive 






r( Q ±+f 6 ± + 3) 


r( «±+^ b ± + 3) 


r( Q ±+f 6 ± + J) 


r( «±+^ b ± + i) 




x tan 1 7T 


a+ + iw6+ 1 

— + - 

2 4 



(13) 



This leads to the exact result 

- — a± —iujb± 



V4(o)_ 
^±(0) ^b ± r( 



2_ r(- 



|)r(2±^i + |) 



-a± —iub± 



x tan it 



2 



- , km- 

iuib± 1 



a± — iub± 



x tan 7r 



2 4 

-ay. + icjfe± 



(14) 



If w has a large positive imaginary part, then the Gamma 
function arguments above tend towards the positive real 
axis, a region where the Gamma function is smooth — 
all potential poles in the logarithmic derivative have been 
isolated in the trigonometric functions. We now use the 
trigonometric identity 



tan A tanB 



to rewrite this as 



cos(A - B) - cos(A + B) 
cos(A - B) + cos(A + B) 



(15) 



V4(Q) 

V±(0) 



2_ r( 



-a± —iujb± 



3\TfO>±-iujb±_ , 3\ 



[)r(- 



cos(7ra±) + sm(iTTtob±) 



a+ —iub^ 



cos(7ra±) — sin(z7rw6±) 
The exact junction condition we wish to apply is 

#h(0) = C(0) 

V>+(o) V-(o)' 



(16) 



(17) 



As long as we are primarily focussed on the highly 
damped QNFs (Im(uj) — > oo) we can employ Stirling's 
approximation to deduce 

I>+|) 



r(z) 

Therefore 

-p/ ±a± —iu:b± 
1 V 2 



1 + 



Re{z) 



(18) 



r( 



±a± —iuib± 



Im(uj)bi 



1 + 



1 



Im(ujb±) 



(19) 



This allows us to deduce an approximate junction condi- 
tion for the asymptotic QNFs 

cos(7ra_|-) + sin(z7ro;6+) cos(7ra_) + sin(i7rw6_) 

cos(7ra + ) — sin(i7ru;6 + ) cos(7ra_) — sin(i7ixjfe_) ' 

(20) 
which is accurate up to fractional corrections of order 
O (1 / Im(u>b±)) . This asymptotic QNF condition can be 
rewritten in any one of the equivalent forms: 

sin(— iirujb+) sin(— iirub-) = cos(7ra + ) cos(-7ra_); (21) 



sinh(7ra;6+) sinh(7ru;6_) = — cos(7ra + ) cos(7ra_); (22) 

cos(— iiruj[b + — &_]) — cos(— i'KLo[b + + &_]) 

= 2 cos(-7ra + )cos(7ra_); (23) 



cosh(7ru;[6 + — &_]) — cosh(7ru;[& + + &_]) 

= 2 cos(7ra + ) cos(7ra_). (24) 

Now suppose 6-f/6_ is rational, that is 

b + /b_=p + /p_, (25) 

and suppose we define 6* by 

b + = p+b*. &_=p_6 >f ; 6* = hcf(& + ,fo_), (26) 

then the asymptotic QNF condition is given by 

sin(— ibjirp+b*) sin(— iujirpJb^) = cos(7ra + ) cos(7ra_). 

(27) 
If W* is any specific solution of this equation, then 

in ( 1 1 \ 

U) n = cj, + — = w* + in 1cm - — , - — (28) 

o» \b + b- ) 

will also be a solution. To characterize all the solutions, 
consider the set of QNFs for which 

Im(uj) < 1/6*, (29) 

and label them as 

uj^ a a e {1,2,3... N}. (30) 

Then the set of QNFs decomposes into a set of "families" 

0J n ,a = ua.a + in/h; (31) 

with a € {1, 2,3... N} and n € {0, 1,2,3.. .}, and where 
N is yet to be detemined. But for rational b+/b- we can 
rewrite the QNF condition as 

cos(— iu>irb*\p+ — P-\) — cos(— iu>nb*\p+ + p~]) 

= 2 cos (7ra+) cos (7ra_) . (32) 



Now define z = exp(w7r6*), then the QNF condition is 

_ z Ip++p-\ _ z -\p++p-] 
4cos (7ra+) cos (7ra_) . (33) 



z \p+-p-\ _|_ z -\p+-p- 



Equivalently 

z 2[p++p-] _ z 2p max + 4cQS ^ a+ ^ cog ( ?ra _- ) Z + [P++P-] 

_ z 2 Pmin +1 = Q (34) 

This is a polynomial of degree N = 2(p + + p_), so it 
has exactly N roots z a (occurring in complex conjugate 
pairs). Thus the QNFs are 



ln(z a ) in 



(35) 



with the imaginary part of the logarithm lying in [0, 27r), 
and where a € {1, 2, 3 ... N} and n e {0, 1, 2,3,.. .}. 

So for rational b+/b- with b+/b- = p+/p- we have 
exactly TV = 2(p + + p~) families of equi-spaced QNFs 



all with with gap i/b* and with (typically distinct) off- 
sets ln(z a )/(7r&»). That is: Arbitrary rational ratios of 
b + /b- automatically imply the uj n = (offset) + in (gap) 
behaviour. 
Now in contrast suppose &+/6- is irrational, that is 

K =hcf(6+,6_) = 0. (36) 

Then the "families" each have only one element 

Lo Q , a o€{l,2,3...oo}. (37) 

That is, there will be no "pattern" in the QNFs, and they 
will not be regularly spaced. (Conversely, if there is a 
"pattern" then b + /b_ is rational.) Stated more formally: 
Suppose we have at least one family of equi-spaced QNFs 
such that 



uj n = luq + inK, 



(38) 



then b+/b- is rational. 

To see this: If we have a family of QNFs of the form 
given in equation (|38[) then we know that Vn > 

+ cos(— iuJon\b + — b—\ + nKir\b+ — 6_|) 
— cos(— iwo7r[& + + b-] + nKir\b + + 6_|) 

= cos(— iujoir\b + — 6_|) — cos(— iuoir[b + + 6-]). (39) 

Let us write this in the form Vn > 

cos(A + nJ) - cos(B + nL) = cos(A) - cos(B), (40) 

and realize that this also implies 

cos(j4+[n + l]J) — cos(B + [n + l]£) = cos(A) -cos(B), 

(41) 
and 

cos(A+[n + 2]J)-cos(B + [n + 2]L) = cos(A) - cos(B). 

(42) 
Now appeal to the trigonometric identity 

cos(A+[n+2]J)+cos(yl+nJ) = 2cos(J) cos(A+[n+l] J), 

(43) 
to deduce 

cos( J) cos(A + [n+l]J)- cos(L) cos(B +[n+ 1]L) 

= cos(A) - cos(B). (44) 



That is, Vn > we have both gj) and fl44|). The first of 
these equations asserts that all the points 



(cos(A+[n + l]J), cos(B + [n+l]L] 



(45) 



lie on the straight line of slope 1 that passes through the 
point (0,cosi3 — cos A). The second of these equations 
asserts that these same points also lie on the straight 
line of slope cos( J)/ cos(L) that passes through the point 
(0, [cos-B — cos A]/ cos L). We then argue as follows: 



|i) If cos J ^ cosL then these two lines are not par- 
allel and so meet only at a single point, let us call it 
(cos A*, cos B*), whence we deduce 

cos(j4+ [n+ I] J) = cos A*; cos(B+ [n+ 1]L) = cosB*. 

(46) 
But then both J and L must be multiples of 2n, and so 
cos J = 1 = cos L contrary to hypothesis. 
|ii) If cos J = cosL ^ 1 then we have both 

cos(A+[n+l]J)-cos(B + [n + l]L) = cos(A) -cos(B), 

(47) 
and 

cos(J) [cos{A + [n + 1] J) - cos(B + [n + 1]L)] = 

cos{A)-cos(B). (48) 

but these are two parallel lines, both of slope 1, that 
never intersect unless cos (J) = 1. Thus cos J = 1 = cos L 
contrary to hypothesis. 

|iii) We therefore conclude that both J and L must be 
multiples of 2-7T, so that in particular cos J = 1 = cos L 
(in which case the QNF condition is certainly satisfied). 
But now 

\b + -b-\/(b + + b_) = J/L G Q, (49) 

and therefore 



b+/b- E Q. 



(50) 



That is: Rational ratios ofb+/b^ are implied by the w„ = 
(offset) +m (gap) behaviour. Thus we have demonstrated 
that the u>„ = (offset) + in (gap) behaviour is generic but 
not universal, and is intimately related to the rationality 
(or otherwise) of the ratio of the e- folding parameters b± . 
Wc have also checked that the analysis sketched above 
satisfies several appropriate consistency checks and has 
suitable well-behaved limits 39]. A particularly impor- 
tant case (not dealt with in 39]) is to consider the situ- 
ation where one side of the potential exhibits power law 
(rather than exponential) falloff. For example, for black 
hole physics a particularly common situation is 



V(x) -> V 0+ 



(x + a) 2 



for x — > +oo. 



(51) 



The exact wavefunctions for this potential can be written 
down in terms of Bessel functions, and in the limit of 
highly damped QNFs one can easily see 



V4(o) 



Im(uj), 



leading to the very simple asymptotic QNF condition 

cos(7ra_) + sin(i7rw6_) 



(52) 



1 = 



cos(7ra_) — sin(z7rw6_) ' 



(53) 



whence, in this particular situation with a one-sided ex- 
ponential falloff one asymptotically has 



(54) 



This observation is useful in that it indicates that one- 
sided exponential falloff can be treated via a minor vari- 
ant of the analysis in [39JJ , and a power law falloff in the 
potential exhibits behaviour qualitatively similar to the 
limit b + — > oo. (As it should on physical grounds.) 

Turning to specific applications in black hole physics: 
Will the general w n — (offset) + in (gap) behaviour dis- 
cussed above extend to more "realistic" astrophysical or 
de Sitter black holes? Consider a "wavepacket" centered 
near the peak of the Regge-Wheeler (Zerelli) potential 
that is built up out of highly damped modes. While the 
initial short-time behaviour of the wavepacket is likely to 
be sensitive to the details of the Regee- Wheeler (Zerelli) 
potential, such a wavepacket will quickly damp out and 
spread out towards both r* — > — oo and r* — > +00, so 
that the wavepacket will penetrate regions where our 
piecewise Eckart model potential should be a good ap- 
proximation to the true potential. We should there- 
fore expect the results of our semi-analytic model to be 
qualitatively (but not necessarily quantitatively) accu- 
rate for estimating the asymptotic QNFs of "realistic" 
black holes. Because of the way the asymptotic QNF 
condition was derived, we do not expect our model to 
give good results for low-lying QNFs. 

Of course one of the key points here is that most of 
the analysis is largely independent of black hole physics, 
and depends only on the falloff conditions placed on the 
potential near asymptotic infinity — indeed if we com- 
pletely forget the black hole motivation, it is already of 
considerable mathematical and physical interest that we 
have a nontrivial extension of the Eckart potential for 
which the QNFs are asymptotically exactly solvable - 
one could in principle loop back to Eckart 's original ar- 
ticle and start asking questions about tunelling proba- 
bilities for electrons encountering such piecewise Eckart 
barriers. 



[1] 
[2] 
[3] 

[4] 



[5] 

[6] 

[7] 
[8] 



C. Eckart, Phys. Rev. 35, 1303 (1930). 
G. Poschl and E. Teller, Z. Phys, 83, 143 (1933). 
P. M. Morse and H. Feshbach, Methods of Theoretical 
Physics, (McGraw-Hill, New York, 1953). 
P. Boonserm, "Rigorous bounds on Transmission, Reflec- 
tion, and Bogoliubov coefficients", PhD Thesis, Victoria 
University of Wellington, 2009, arXiv: 0907.0045 [math- 
ph]. (See especially chapter 3 for a discussion of the con- 
voluted history of the Eckart potential.) 
S. Chandrasekhar and S. Detweiler, Proc. Roy. Soc. 
Lond. A 344 (1975) 441. 

K. D. Kokkotas and B. G. Schmidt, Living Rev. Rel. 2 
(1999) 2 arXiv:gr-qc/9909058 . 
H. P. Nollert, Class. Quant. Grav. 16 (1999) R159. 
E. Berti, V. Cardoso and A. O. Starinets, Class. Quant. 
Grav. 26 (2009) 163001 arXiv:090 O975l [gr-qc]]. 



[9 

[10 

[11 

[12 
[13 

[14; 
[15; 

[is; 

[17 
[18 

[19 

[20 

[21 

[22 
[23 
[24 
[25 

[26 

[27 

[28 

[29 

[30 

[31 

[32 
[33 
[34 
[35 
[36 
[37 
[38 
[39 
[40 
[41 



O. Dreyer, Phys. Rev. Lett. 90 (2003) 081301 

[arXiv:gr-qc/0211076] . 

J. Natario and R. Schiappa, Adv. Theor. Math. Phys. 8 

(2004) 1001 [arXiv:hep-th/0411267| . 

V. Ferrari and B. Mashhoon, Phys. Rev. D 30 (1984) 

295. 

S. Iyer and C. M. Will, Phys. Rev. D 35 (1987) 3621. 

C. M. Will and J. W. Guinn, Phys. Rev. A 37 (1988) 
3674. 

J. W. Guinn, C. M. Will, Y. Kojima and B. F. Schutz, 

Class. Quant. Grav. 7 (1990) L47. 

R. A. Konoplya, Phys. Rev. D 68 (2003) 024018 

arXiv:gr-qc/0303052; . 

A. J. M. Medved, D. Martin and M. Visser, Class. Quant. 
Grav. 21 (2004) 1393 arXiv:gr-qc/0310009 . 
A. J. M. Medved, D. Martin and M. Visser, C lass. Quant. 
Grav. 21 (2004) 2393 arXiv:gr-qc/0310097 . 
T. Padmanabhan, Class. Quant. Grav. 21 (2004) LI 
[arXiv:gr-q c/031002~7]. 

T. R. Choudhury and T. Padmanabhan, Phys. Rev. D 
69 (2004) 064033 [arXiv:gr-qc/0311064| . 
N. Andersson, Proc. Roy. Soc. Lond. A 439 (1992) 47- 
58. 

N. Andersson and C. J. Howls, Class. Quant. Grav. 21 
(2004) 1623 arXiv:gr-qc/0307020; . 

E. W. Leaver, Proc. Roy. Soc. Lond. A 402 (1985) 285. 
E. W. Leaver, Phys. Rev. D 34, 384 (1986). 
E. W. Leaver, Phys. Rev. D 41 (1990) 2986. 
L. Motl, Adv. Theor. Math. Phys. 6 (2003) 1135 

arXiv:gr-qc/0212096 . 
L. Motl an d A. Neitzke, Adv. Theor. Math. Phys. 7 

(2003) 307 arXiv:hep-th/0301173;. 

S. Das and S. Shankaranarayanan, Class. Quant. Grav. 
22 (2005) L7 arXiv:hep-th/0410209; . 

A. Ghosh, S. Shankaranarayanan and S. Das, Class. 
Quant. Grav. 23 (2006) 1851 arXi v:hep-th/0510186| . 

J. S. F. Chan and R. B. Mann, Phys. Rev. D 55 (1997) 

7546 arXiv:gr-qc/9612026| . 

J. S. F. Chan and R. B. Mann, Phys. Rev. D 59 (1999) 

064025. 

G. T. Horowit z and V. E. Hu beny, Phys. Rev. D 62 

(2000) 024027 arXiv:hep-th/9909056 . 

B. Wang, C. Y. Lin and E. Abda lla, Phys. Lett. B 481 

(2000) 79 arXiv:hep-th/0003295; . 

B. Wang, C. M olina andTTH bdalla, Phys. Rev. D 63 

(2001) 084001 [arXiv:hep -th/0005143 . 

J. M. Zhu, B. Wang and E. Abdalla, Phys. Rev. D 63 

(2001) 124004 arXiv:hep-th/0101133 . 

V. Cardoso and J. P. S. Lemo s, Phys. Rev. D 64 (2001) 

084017 [arXiv:gr-qc/0105103j . 

B. Wang, E. Abdalla and R. B. Mann, Phys. Rev. D 65 

(2002) 084006 [arXiv:hep- th/0107243 . 

D. P. D u, B. Wang and R. K. S u, Phys. Rev. D 70 (2004) 
064024 arXiv:hep-th/0404047 . 

B. Wang, C. Y. Lin and C. M olina, Phys. Rev. D 70 

(2004) 064025 [arXiv:hep -th/0407024 . 

J. Skakala and M. Visser, "Semi-analytic results for 
quasi-normal frequencies," arXiv:1004.2539 [gr-qc]. 
V. Suneeta , Phys. Rev. D 68 (2003) 024020 
arXiv:gr-qc/0303114 . 

H. R. Beyer, C ommun. Math. Phys. 204 (1999) 397 
[arXiv:gr-qc/9803034j . 



